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Within the minimal supergravity model (mSUGRA) framework, the expectation for the relic 
density of neutralinos exceeds the WMAP determination, unless neutralinos a) have a significant 
higgsino component, 6) have a mass close to half that of a heavy Higgs boson, or c) can efficiently 
co-annihilate with a charged or colored particle. Within a 1-parameter extension of the mSUGRA 
model which includes non-universal Higgs masses, we show that agreement with the WMAP data 
can be obtained over a wide range of mSUGRA parameters for scenarios a) and 6), so that the 
phenomenological implications may be much more diverse than in mSUGRA. We show that direct 
and/or indirect detection of neutralino dark matter should be possible at various current and planned 
facilities. 
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Supersymmetry (SUSY) is a novel spacetinie symme- 
try between bosons and fermions. In realistic models, 
SUSY is broken at the weak scale implying that all Stan- 
dard Model (SM) particles must have superpartners with 
masses in the range ^ 100 — 1000 GcV that will be acces- 
sible to colliders. The lightest neutralino Zi of i?-parity 
conserving SUSY models is an especially attractive and 
well-motivated candidate for cold dark matter (CDM) in 
the umverse[jy|. The WMAP collaboration,^ has recently 
determined the relic density of CDM in the universe to 
be QcDAih"^ = 0.113 ±0.009 (Icr); this measurement, es- 
pecially its implied upper limit, provides a powerful con- 
straint for any model of particle physics that includes a 
candidate for CDM. 

The present neutralino relic density can be deter- 
mined by solving the Boltzmann equation for neutrali- 
nos in a Friedmann-Robertson- Walker universe. The 
central aspect of the calculation involves a compu- 
tation of the thermally averaged neutralino annihila- 
tion and co-annihilation cross sections. Many analy- 
ses of SUSY CDM have been carried out within the 
paradigm minimal supergravity (mSUGRA) modelQ- 
The free parameters of the mSUGRA model consist of 
mo, mi/2, Aq, tan/9, and s«(7n(/i), where mo is the com- 
mon scalar mass, mi/2 is the common gaugino mass, and 
Aq is a common trilinear soft SUSY breaking parameter 
all defined at the scale Q = Mqut 10^^ GeV. The 

parameter tan (3 is the ratio of weak scale Higgs field vac- 
uum expectation values, and is a superpotential param- 
eter whose magnitude is constrained by the requirement 
of radiative electroweak symmetry breaking (REWSB). 
Once these model parameters are specified, then all spar- 
ticle masses and mixings are determined, and scattering 
cross sections, decay rates, relic density and dark matter 
detection rates may all be reliably calculated. 

It has been increasingly recognized that the typical 
value of Q~ h'^ as given by the mSUGRA model signifi- 



cantly exceeds its WMAP upper limit. Only specific re- 
gions of the mSUGRA model parameter space where neu- 
tralinos can efficiently annihilate are in accord with the 
WMAP data. These include: i) The bulk annihilation re- 
gion at low values of mg and »tii/2, where neutralino pair 
annihilation occurs at a large rate via i-channel slepton 
exchange, ii) The stau (stop) co-annihilation region at 
low Too 4] (for special values of vlo'sj) where to— ~ ruf-^ 

(jn~ ~ Tif^) so that Zis may co-annihilate with fis (ti) 
in the early universe. Hi) The hyperbolic branch/focus 
point (HB/FP) region at large Too '^3 — 7 TeV (de- 
pending sensitively on the assumed value of mt) near the 
boundary of the REWSB excluded region where be- 
comes small, and the neutralinos have a significant hig- 
gsino component, which facilitates annihilations to WW 
and ZZ pairs j^]. iv). The vl-annihilation funnel, which 
occurs at very large tan/3 ~ 45 — 60. In this case, the 
value of niA ~ 2to~ , so that neutralino annihilation in 
the early universe is enhanced by the A (and also H) 
s-channel poles 0. These allowed regions frequently oc- 
cur either at edges of parameter space, or for extreme 
values of mSUGRA parameters. In this paper, we show 
that mSUGRA parameter choices that were previously 
disallowed by WMAP can in fact be brought into accord 
with the relic density measurement via two possible so- 
lutions, in a well motivated one parameter extension of 
the model. Thus, parameter choices thought to be irrele- 
vant for collider and other SUSY searches are in fact now 
allowed. 

The assumption of a universal soft SUSY breaking 
scalar mass too is phenomenologically motivated for mat- 
ter scalars, since universality guarantees in a super-GIM 
mechanism which suppresses dangerous flavor chang- 
ing neutral current processes Q. Within the frame- 
work of gravity-mediated SUSY breaking, universality 
for all scalars is realized by assuming a flat Kahler po- 
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tential, though non-universal scalar masses are certainly 
possible [^|. While universal masses for matter scalars 
is phenomenologically motivated, there is, no motivation 
for assuming that the Higgs field scalars also have a com- 
mon GUT scale mass rriQ. Indeed, in supersymmetric 
SO{10) grand unified theories, which are highly moti- 
vated by recent data on neutrino masses, the matter mul- 
tiplets of each generation belong to the 16 dimensional 
spinorial representation of 50(10), while the Higgs mul- 
tiplets inhabit a single 10 dimensional fundamental repre- 
sentation, and would have an unrelated mass, assuming 
that the SUSY breaking mechanism is not completely 
blind to gauge quantum numbers. In the simplest case, 
we may naively expect the Higgs soft SUSY breaking 
squared masses to obey m^^^ — ^ niQ so that this 
one parameter extension is characterized by, 

mo, m^, mi/2, ^o, tan/?, and sign{^), (1) 



m„=1000GeV, m„^=200GeV, tanP=20, A„=0, |i>0, m=178GeV 
m„=300GeV, m|,2=300GeV, tanp=10, A„=0, |J,>0, m=178GeV 



where 



We dub this the 



non-universal Higgs mass (NUHMl) model, to distin- 
guish it from the more general NUHM2 model0, 0| 
where each Minimal Supersymmetric Standard Model 
(MSSM) Higgs scalar has an independent mass param- 
eter at Q = Mgut- The NUHM2 model has recently 
been investigated by Ellis et al. and is motivated in- 
stead by SU{h) SUSY grand unification. These NUHM2 
models have also been investigated recently within the 
context of Yukawa unified SUSY models 0. 

We begin our analysis of the NUHMl model by not- 
ing that the parameter range for need not be limited 
to positive values. Some authors have invoked the so- 
called GUT stability bound, sm'^^ + ^'^{Mgut) > (with 
s = siqn{m^)), to avoid EWSB at too high a scale; see 
Ref. IUJI where the applicability of these bounds is dis- 
cussed. Except to note that these are hardly constrain- 
ing for the cases that we present here, we do not consider 
them any further. 

We generate sparticle mass spectra using ISAJET 
v7.72[l^ upgraded to allow the input of negative Higgs 
squared masses. ISAJET includes full one-loop radiative 
corrections to all sparticle masses and Yukawa couplings, 
and minimizes the scalar potential using the renormal- 
ization group improved effective potential renormalized 
at an optimized scale choice to account for leading two 
loop terms. We evaluate the relic density using the 



^ We stress that we are using 50(10) considerations only to guide 
our thinking, and do not commit to any particular SO(IO) frame- 
work. As a result, many of the usual consequences such as 
the unification of Yukawa couplings that follow in the minimal 
SO(IO) model, no longer obtain in our case. In particular, as 
we will see below, it is possible to break electroweak symmetry 
radiatively even if the MSSM Higgs scalar doublets have equal 
mass parameters at Q = Mqux- 
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FIG. 1: The values of a), fi, tua and m~ versus m^t/mo for 

two cases in the NUHMl model. We also show h). Q,cDMh^ , 
c). BF{b — > 57) and d). Aa^. The curves terminate on the 
left (right) because m\ < (p^ < 0) so that REWSB is not 
obtained. 



IsaReD[T^ program, and to evaluate the indirect CDM 
signals expected from the NUHMl model, we adopt the 
DarkSUSYJJi] package interfaced to ISAJET. 

As a first example, in Fig. Qwe show one of the crit- 
ical aspects of the NUHMl model, where we plot in a) 
the values of /i, tua and m~ versus m^/wo, while fix- 
ing other parameters as listed. One can see that fj, be- 
comes much larger for < —mo, and much smaller for 
TO0 > mo compared to the case of m^ = mo. The region 
of small fi is of particular interest since in that case the 
lightest neutralino becomes more higgsino-like, and gives 
rise to a relic density which can be in accord with the 
WMAP determination [1^. In the mSUGRA model, the 
higgsino-like Zi region occurs in the HB /FP region men- 
tioned above, with mo ~ 3 — 7 TeV. We immediately see 
one important virtue of the NUHMl model: the higgsino 
annihilation region may be reached even with relatively 
low values of mo and m^. Thus, unlike the mSUGRA 
model case, the low region of the NUHMl model can 
occur for small values of mo , and mx/21 and so need 
not suffer from fine tuning[l9j. 

We also see from Fig. ^ that the value of tua can 
range beyond its mSUGRA value for large values of m^, 
to quite small values when m^ becomes less than zero|20j. 
In particular, when tua ~ 2m~ , then neutralinos in the 
early universe may annihilate efficiently through the A 
and H Higgs resonances, so that again flcDMh^ may be 
brought into accord with the WMAP determination. In 
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the mSUGRA model, the ^-annihilation funnel occurs 
only at large tan/3 > 45 (55) for ^ < (^ > 0). How- 
ever, in the NUHMl model, the A-funnel region may be 
reached even for low tan/3 values, if is taken to have 
negative values. At low tan/3, the A and H widths are 
much narrower than the large tan /? case. We also show 
in frames b) the relic density fi/i^, c) the branching frac- 
tion BF{b sj) (allowed range conservatively taken to 
be - (2.6-4.5) x 10"^), and d) the SUSY contribution to 
the muon magnetic moment Aa^. We have also checked 
that the branching ratio B{Bs fi^n^) < 10~* even 
for the smallest values of tua in Fig. ^ 

To understand the behavior of the /i parameter 
and niA in the NUHMl model, we define Amj^ ^ = 
mfj^ ^ (NUHMl )-r7i|j^ ^(mSUGRA), and obtain the one 
loop rcnormalization group equation 



dAm'jj 
Jt 



X 3/2 [Xt (NUHMl) -Xt(mSUGRA)] 



(2) 



where t = log(Q), ft is the top quark Yukawa coupling, 



and Xt — 
integrated to obtain. 



■ At . Eq. 10) can be readily 



Amjj {weak) 



Amjj (GUT) X e""'* , 



where Jt — J dtft > 0. We see that Amjj maintains 
its sign under KG evolution, and furthermore, reduces in 
magnitude in the evolution from high to low scales. The 
same considerations also apply to Amj^^ , except that the 
effect is much smaller because for the modest values of 
tan/3 that we consider, ft^r ^ ft so that Jb and Jr are 
both < 1 and Am'\j^{GUT) ~ Amj^^iweak). From the 
tree level minimization condition for EWSB in the MSSM 



^ P7 



tan^/? 



(tan^ /? - 1) 



we see that for moderate to large values of tan/3 (as 
favored by LEP2 Higgs boson mass constraints), and 

I^'hJ ^ ^^i' ^"^Ifu- ^^us t^^t ^ large 

value of fj? is expected if is large and negative. In 
contrast, for > mg, we expect a smaller weak scale 
magnitude of mj^ , and hence of /i^ . 

The tree level pseudoscalar Higgs mass ruA is given by 



'Ha 



where the last equality holds in the approximation /i^ 



'H^ 



For values of tan/3 such that Jb and Jr can be 



neglected compared to Jt, it is easy to see that 



TO^ (NUHMl)- 



(mSUGRA) = (sml-mDil-e- 



thereby accounting for the uia dependence in Fig. 

Most other sparticle masses are relatively invariant to 
changes in (justifying our approximation in the last 



step of 0). An obvious exception occurs for the chargino 
Wi and neutralino Zi^2 masses, which become small when 
fi ~ Ml 2, and the Zi becomes increasingly higgsino-like. 
The other exception occurs for third generation squark 
masses. In this case, the (^Li^l) and in running masses 
depend on to^ via the ft^t terms in the RG equations, 
and can be somewhat suppressed for large positive values 
of ml. 

In order to estimate the dark matter detection rates 
for the NUHMl model, we adopt a cored halo model 
(the Burkert profile, j21|) which has been tested against 
a large sample of rotation curves in spiral galaxies |22|. 
and whose density profile reads 



PB{r) 



(l + r/a)(l + (r/a)2)' 



(3) 



with a = 11.7 kpc and p% = 0.34 GeV cm"^. The par- 
ticular configuration we use has been found after imple- 
menting all available dynamical constraints and numeri- 
cal simulation indications on the halo mass-concentration 
correlation |23!|. The corresponding velocity distribution 
has been self-consistently computed (see Ref. [23] for de- 
tails), thus allowing for a reliable comparison among di- 
rect and indirect detection techniques ^24]. 

Following Ref. ,2^], we adopt here Visibility Ra- 
tios (VR), i.e. ratios of the expected signals from a 
given supersymmetric model over the projected future 
sensitivities in the particular detection technique, and 
for the corresponding neutralino mass. Direct detec- 
tion refers to the spin-independent (SI) neutralino-proton 
scatterin g cr oss section, compared to the sensitivity of 
CDMS-Ill2a| and of future ton-size experiments, such as 
XENON [27|. We compute the flux of neutrino- induced 
muons produced by ncutralino pair annihilations in the 
sun, and compare the flux against the projected sensi- 
tivity of the km^-size detector IceCube (computed for 
an energy threshold of 1 GeV), taking into account the 
energy threshold mismatch and the dependence of the 
detector sensitivity on the soft {e.g. bb, in the left 
panel of Fig. 01 and hard {e.g. W'^W^ , in the right 
panel of Fig. |2Jl neutrino production channels, follow- 
ing Ref. 28] . As regards antideuterons, we compute the 
solar-modulated flux, at solar maximum, in the low en- 
ergy range 0.1 < Tjj < 0.4 GeV, where the background 
is strongly suppressed ^2^. The VR is then defined as 
the ratio of the computed D flux over the sensitivity of a 
GAPS ^3 detector placed on a satellite orbiting around 
the Earth, and tuned to look for antideuterons in the 
mentioned very low kinetic energy interval, correspond- 



ing to a D flux of 2.6 x lO^-'m-^sr^^GeV^'s" 



r resp c 

13. 



Re- 



garding antiprotons and positrons, we adopt the quantity 
/0, defined in Ref. [3], and examine the sensitivity of the 
PAMELA experiment after 3 years of data-taking. Fi- 
nally, the gamma ray flux from the galactic center (GC) 
refers to the sensitivity of the GLAST experiment to the 
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tan|3=10, |J.>0, Aj,=0, 01^=300 GeV 
m,„=3()0 GeV, m =178 GeV 
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FIG. 2: Dark matter detection Visibility Ratios {i.e. signal-to-projected-sensitivity ratios) for various direct and indirect 
techniques, as a function of the GUT-scale non-universal Higgs mass parameter m^. A Visibility Ratio larger than 1 implies 
that the model will be detectable. The dark (red) shadings to the left (right) of the plot indicate regions ruled out by the LEP2 
bounds on tua ("^^ )• Within the lighter shaded (green) regions a thermal neutralino relic abundance compatible with the 



WMAP upper bound ^ is produced. 



integral photon flux above a 1 GeV threshold |3ll |. 

While only large values of will be probed via CDM 
detection at CDMS-II, all the parameter space will be 
probed by future experiments such as XENON. Inter- 
estingly enough, we obtain observable rates for indirect 
detection experiments only in those regions where a neu- 
tralino relic abundance compatible with the WMAP up- 
per bound is produced. The first cosmologically al- 
lowed parameter space region, lying in the range —2.57 < 
m^/rriQ < —2.45, features a large neutralino annihila- 
tion cross section, due to the proximity of heavy Higgs 
resonances 3^ . This yields detectable rates in all exper- 
iments apart from IceCube, where the rate of neutralino 
capture inside the sun depends critically on the spin- 
dependent neutralino-nucleon coupling, which is sup- 
pressed in this large region: see Fig. |2 In the second 
viable parameter space region (1.83 < m^/mo < 1.93), 
a large higgsino fraction implies instead detectability 
at IceCube, as well as at future antiprotons and an- 
tideuterons searches. Though positrons and gamma rays 
rates are enhanced in this region, they will still be slightly 
below the respective planned experimental sensitivities. 
We note, however, that detection rates from halo annihi- 
lations are very sensitive to the CDM halo profile which 
is assumed. More quantitatively, we find that with an 
extremely cuspy profile (such as the adiabatically con- 
tracted NFW profile of Ref. H^), the direct detection and 
neutrino telescope rates increase only by around 20%, 
the D flux by ^ 10^, positrons and antiprotons VR's by 
^ 10^, and the gamma ray flux from the GC by more 
than 10^ (see also the discussions in Ref. [s^ and [2^). 

Summary: The mSUGRA model prediction of the 
CDM relic density can be brought into accord with the 
WMAP measured value only if the neutralino annihila- 
tion cross section is enhanced, usually due to A and H 
boson resonances, or by increasing the higgsino content of 



the neutralino. Within the mSUGRA framework, these 
possibilities occur only if tan/3 is very large or scalar 
masses are in the HB/FP region at very large toq. We 
have shown that in the NUHMl model, a well motivated 
1-parameter extension of mSUGRA, neutralino annihila- 
tion via heavy Higgs resonances can occur even for low 
values of tan/?, and neutralino annihilation via higgsino 
components can occur at low values of tuq. Thus, a 
much wider range of collider signals and low energy phe- 
nomenology is possible compared to expectations from 
the mSUGRA model. 
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